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Abstract 

The Hamiltonian describing a single ion placed in a potential trap in in- 
teraction with a laser beam is studied by means of a suitable perturbative 
approach. It is shown, in particular, that the rotating wave approximation 
does not provide the correct expression, already at the first perturbative 
order, of the evolution operator of the system. 



1 Introduction 

Trapped ions in interaction with laser beams are extremely useful tools for inves- 
tigating fundamental aspects of quantum physics. For instance, they have been 
used for the generation of coherent, squeezed PP and Schrodinger-cat states 0, 
and for the preparation of entangled Bell and GHZ states [3] 0]. They have 
also had important experimental applications as, for instance, precision spec- 
troscopy [5] and laser cooling El- 

Recently, the interest for laser-driven ion traps has received a novel impulse 
in view of its applications in the fastly developing area of quantum comput- 
ing. Indeed, in a quantum computer (QC), information is stored in a 'quantum 
register' composed of N two-level systems representing the quantum bits, or 
qubits The storage of data and all the basic operations are implemented 
by inducing controlled dynamics on the quantum register |1(J) . Since a QC is 
a quantum mechanical system, it can perform superpositions of computation 
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operations with a remarkable gain of efficiency with respect to a classical com- 
puter; a typical example is the solution of the problem, fundamental for some 
cryptographic schemes, of factoring large numbers into primes In 1995 

Cirac and Zoller ^21 proposed a concrete model for a ion-trap computer con- 
sisting of N atomic ions trapped in a parabolic potential well, each ion being 
regarded as a two-level system, hence as a realization of a qubit. The control of 
the quantum degrees of freedom is achieved by addressing the ions with time, 
frequency and intensity controlled laser beams (see also ^31 El)- 

All the applications mentioned above are realized by ion traps with a linear 
geometry. In such traps a strong confinement is induced along the y and z axes, 
while a weak harmonic confinement is induced along the principal trap axis x. 
This confinement scheme is realized by the Paul linear trap ^BJ. N ions in a 
linear trap will form a chain along the principal trap axis. This allows to reduce 
the initial 3iV-dimensional model to a A-dimcnsional one, which can be treated 
conveniently by the introduction of the normal coordinates of the ion chain [16| . 
As already mentioned, controlled dynamics can be induced on the ion trap by 
means of laser beams. The Hamiltonian describing the total system \IQ \H\ ~ 
which we will call simply the 'ion trap Hamiltonian' (ITH) — in spite of its rela- 
tive formal simplicity, gives rise to a Schrodingcr equation whose exact solutions 
are not known and its study requires the adoption of suitable approximations. 
The stantardly used method is the rotating wave approximation (RWA) 
|19|. The RWA is a very popular technique in quantum optics and, in general, 
in the study of resonance phenomena since it leads to considerable simplifica- 
tions in many calculation procedures. It consists essentially in passing to the 
interaction picture and then dropping those terms of the effective Hamiltonian 
which are rapidly oscillating (usually called 'counter rotating terms' or 'virtual 
terms'). Using this technique and the Lamb-Dicke approximation [201, the ITH 
can be reduced to an effective Hamiltonian formally identical to the Hamilto- 
nian of the Jaynes-Cummings model (JCM) |2J EH in the interaction picture, 
hence explicitly integrable. 

Despite its great popularity, the general validity of the RWA is rather uncer- 
tain. In particular, the application of the RWA to the evaluation of the evolution 
operator is justified only by semiqualitative arguments and the ubiquity of this 
procedure is mainly due to the chance of performing explicit calculations par- 
tially supported by the prediction of some experimentally observable phenom- 
ena, for instance the typical 'collapses and revivals' in two-level systems (23] HH 
|25| . Thus, some attempts of taking into account the impact of counter rotating 
terms have been made. Specifically, perturbative corrections to the energy spec- 
trum and corrections to the time evolution by means of path integral [2Zj 
and perturbative |2E] QUI EH] techniques have been investigated. It should be 
also mentioned that an attempt of considering the counter rotating contribu- 
tions is already present in the classical study of the magnetic resonance done by 
Bloch and Siegert |3J and in the later related work of Shirley |22] ■ Anyway, the 
validity of all these approaches, as well as of the RWA itself OH, rests on the 
smallness of the coupling constant, which in the case of a laser-driven ion trap 
is proportional to the intensity of the laser field. This is a severe drawback since 
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an intense laser field implies a fast coupling between the two internal energy 
levels of the trapped ions i.e., for instance, a fast QC. 

In the present paper, we propose to study the ITH using a perturbative ap- 
proach along new lines with respect to the exhisting literature. Since, by virtue 
of the normal coordinates of the ion chain, the iV-ion case does not introduce 
any essential complication with respect to the single ion case (especially if the 
nonharmonic component of the ion- ion interaction can be neglected) , we will re- 
strict, for the sake of simplicity, to the latter case. Obviously, the TV-ion case is 
of great interest for the applications and will deserve a particular attention in a 
forthcoming paper We proceed as follows. First, we notice that passing to 
a 'rotating frame', i.e. to a suitable interaction picture, the ITH transforms into 
a time-independent Hamiltonian which we will call the 'rotating frame Hamilto- 
nian' (RFH). At this point, we use the fact that — as it has been shown by some 
of the authors [jSj — the RFH is unitarily equivalent to a Hamiltonian formally 
similar to the RFH, except for the fact that the new coupling constant is not 
proportional any more to the field intensity but a simple bounded function of it 
(see section|2}. We will call this Hamiltonian the 'balanced Hamiltonian' (BH). 
The BH is an ideal starting point for a perturbative approach, since the results 
obtained by its study hold also in the strong field regime. Next, using the tools 
of perturbation theory for linear operators |37| |38|. we develop a perturbative 
procedure which allows to give approximate expressions of the evolution opera- 
tor associated with the BH (hence with the ITH) in terms of unitary operators. 
Our basic idea is the following. The Jaynes-Cummings Hamiltonian is exactly 
solvable — i.e. its spectrum is discrete and its eigenvalues and eigenvectors are 
known — by virtue of the fact that it can be written as the sum of two com- 
muting operators and one of the two is trivially solvable (see section EJl. Then, 
in general, we wonder if, given a perturbed Hamiltonian = -fio + Ajjt, with 
$3o exactly solvable, we can build a procedure which allows to compute, for each 
n = 1,2, .. ., hermitian operators ijg (A), £^(A) depending analytically on the 
parameter A, such that 

53(A) =^ n) (A) + c(")(A) + 0(A" +1 ), (1) 

5}q TI ' ) (A) is exactly solvable and 

-0. (2) 

Indeed, in this case the determination of an approximate expression, at each 
perturbative order, of the evolution operator associated with .f)(A) is greatly 
simplified. We will show in section 0] that this idea is correct and we will 
apply the method to the BH in section |SJ In sections and [7| we compare the 
perturbative expressions obtained for the energy spectrum and the evolution 
operator with the results obtained applying the RWA to the BH. It will be shown 
that, while the RWA gives the right first order correction to the unperturbed 
eigenvalues, it does not give the right first order expression for the evolution 
operator since it neglects corrections to the eigenprojectors. Eventually, in 
section 03 conclusions are drawn. 



£(")(A),^ n) (A) 
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2 The Jaynes-Cummings model 



In this section we will give a very concise treatment of the classical JCM. This 
will allow us to gain a better insight in our perturbative analysis of the ion trap 
Hamiltonian. 

The Hamiltonian of the JCM reads (h = 1): 

Hj C = cn+-w(T z +Ai/(a(j + +a t (T_), (3) 

where a, , n are the annihilation, creation and number operators and a z , o± 
the Pauli operators. In the following, we will denote by {\n}} the Fock basis 
and by \g) , |e) the eigenvectors of a z : 

h\n)=n \n) , a z \g) = - \g) , <r z \e) = |e) . 

The JCM is exactly solvable. This is due to the fact that its Hamiltonian can be 
easily represented as the sum of two constants of the motion. Indcccd, observe 
that Hjc can be written as Hjc = Af + S, where 

Af := v {n + , S := ^ (u) - v) o z + \v (a<r + + a f er_) (4) 

and 

[Af,S] = [Af,H JC } = [S,Hjc] = 0. (5) 
It follows that the evolution operator of the JCM factorizes as 

e -iH Jc t _ e -iNt e ~iSt 

Moreover, Af has a discrete spectrum and its eigenspaces, namely the one- 
dimensional eigenspace Ho = Span{|0)<g)|g)} and the two-dimensional eigenspaces 

Tin — Span{|n - 1) ® |e), |n) ® \g)}, n=l,2, ... , 

are invariant subspaces for S, which can be diagonalized in each of these mutu- 
ally orthogonal subspaces. This yelds to a diagonalization of the whole Hamil- 
tonian and to an explicit expression for the unitary operator e~ tSt . If the 
resonance condition v — u is satisfied, then this operator assumes a particularly 
simple form and can be also computed by direct exponentiation, indeed, using 
the fact that 

{ao- + +a'<7-) = [a a' a + er_ + a' a er_ cr + J 

= h m \g)(g\ + (h+ir\e)(e\, 

we find: 

JJ(t) := exp(-iSt) 

= cos v Vft tj \g) (g\ + cos v Vn + 1 tj 

, sin (A v VnTT t) sin ( A v *) . 
i V - = '- ao-++ '- ci (T_ I . (6) 

V n + 1 Vn 



e e 
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We will call JO(t) the Jaynes-Cummings evolutor. 
We conclude this section observing that since 

e^"«e'^" = ia, 

where e 1 ? " is nothing but the Fourier-Planchcrel operator, the Hamiltonian iJjc 
is unitarily equivalent to the following one: 

#jc := e"^ ft Hj C e^ h =N + S, (7) 

where ^ 

5 := - (uj — v) a z + i A v [a a + — a) a J) . (8) 



3 The ion trap Hamiltonian 

A two-level ion of mass \x in a potential trap, with strong confinement along the 
y and z axes, and weak harmonic binding of frequency v along the x-axis (the 
'trap axis'), can be described neglecting the motion of the ions transverse to 
the trap axis — by a Hamiltonian of the following type (h= 1): 



+ 1 
H Q = va>a + -uJ ge (. 

where a is the vibrational annihilation operator 



P., 



and a z the effective spin operator associated with the internal degrees of freedom 
of the ion. Let us suppose now that the ion is addressed by a laser beam of fre- 
quency ojl in a traveling wave configuration. Then, the Hamiltonian describing 
the physical system (ITH) becomes 

H(t) = H + H t (t), (9) 

where 

H t (t) := n R (e^< a- D(irtf + e"^' a+ D(i V )) , (10) 

with £Ir — p£ the Rabi frequency and £ the intensity of the laser field. More- 
over, we have set: 

D(irj) := cxp(ir? (a + a))) , (11) 

where 

kreosd) 

* == W (12) 
- with ki the wavevector and <j> the angle between the a;-axis and Icl — is 
the Lamb-Dicke factor. Notice that D(a), a e C, is a displacement operator, 
namely 

D(a) := exp (a a) — a* a) , D{a) a D{a)^ = a — a. (13) 
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In order to work with operator matrices, the Hilbert space of the total system 
('pointlike' ion + internal degrees of freedom of the ion), namely Tt (g> C 2 , Ti = 
L 2 (R), will be identified with the space U®U. 

Now, the dynamical problem associated with the time-dependent Hamilto- 
nian H(t) can be turned into a time-independent problem. Indeed, switching 
to the interaction picture with reference Hamiltonian ^LULta z and setting 

Rt •= exp (iijULtcr^j , (14) 

one obtains the time-independent 'rotating frame Hamiltonian' (RFH) 

H = R t (H{t)- l -Lo L a^R\ 

= vn+^S(T z +n R (a-D{if])' i +<j+D(ir])) , (15) 

where n — a^a is the number operator and S := ui ge — u>l is the ion-laser detun- 
ing. 

At this point, in order to give an approximate expression of the evolution opera- 
tor associated with the Hamiltonian H , the rotating wave approximation (RWA) 
is usually applied. It consists in expanding the exponential D{irf), then pass- 
ing to the interaction picture with reference Hamiltonian Hq — v n + ^8 o~ z , so 
obtaining the interaction picture Hamiltonian 

H int (t) = n R (l - t V (e wt a f + e~ wt a) + . . .) e~ l5t a- + h.c. , 

and, finally, retaining only that terms in iJj nt (i) which are slowly rotating. 
Hence — assuming that rj <C 1 (Lamb-Dicke regime), so that one can keep only 
the terms which are at most linear in 77 — in correspondence to the three types 
of resonance 

S = uj ge - lo l ~ 0, (5 + i/~0, £-i/~0, (16) 

one obtains respectively the following three types of effective interaction picture 
Hamiltonian: 

fijg = n R (a. +*+), 
= it] £l R (a f a + - a cr_) , H^' = iij Q R (a a + - a f a J) . 

These effective Hamiltonians, in correspondence to the respective resonances, 
commute with the reference Hamiltonian Hq. This is due to the fact that the 
resonances (|16|) are associated with degeneracies of the reference Hamiltonian. 
In fact, it turns out that the spectrum of Hq is degenerate if and only if the 
condition mv = \5\, m = 0, 1, 2, . . ., holds. Notice that, in particular, H^ is 
equal — up to a unitary transformation and setting 5 = u> and 77 £l R — X v — 
to the constant of the motion S of the Jaynes-Cummings Hamiltonian in the 
resonant regime (see section |2J. 
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It will be shown that the RWA is a rather poor approximation for the eval- 
uation of the evolution operator. The argument usually adopted in order to 
support its validity is the following. Let us consider the Feynman-Dyson expan- 
sion of the evolution operator Ui n t(t, to) associated with the interaction picture 
Hamiltonian Hi nt (t). At the first order one has: 



U iat (t, t ) ^ Id -i f flint (t) dr. 

Jta 



to 

Then, it is argued that the fastly oscillating terms in ffi n t give a smaller contri- 
bution to the integral in the r.h.s. of the previous formula with respect to the 
slowly rotating ones. As it will be seen later on, this argument turns out to be 
erroneous. The misunderstanding stems from the fact that one is using a per- 
turbative expansion of the evolution operator whose terms (except the identity) 
are not unitary operators. 

Besides, any perturbative approach does not work in the strong field regime, 
since the Rabi frequency FIr which appears in the interaction component of the 
RFH is proportional to the intensity of the laser field. This problem can be 
bypassed by means of a suitable unitary transformation of the Hamiltonian H 
(see [SI]). This transformation allows to obtain a 'balanced Hamiltonian' (BH) 
H which is the sum of a large component having a simple 'diagonal' form — 
i.e. such that its matrix representation in the standard basis 

{\n)®\g), |n)o|e) : n = 0,1,2,...} (17) 

is diagonal — and a small interaction component scarcely sensitive to the Rabi 
frequency Or. Indeed, introduced the dimensionless parameter A := 6/SIr 
(from this point onwards we will assume that flji ^ 0), there is a unitary 
operator Ta such that 

H :=T A HTl = H + H h (18) 



where we have set: 



with 



Hq := vh + -^S(T Z + Xfjv, (19) 



:ri, A:= —L= r) = -^rj, (20) 



V4 + A 2 V4 + A 2 A 



and 



#j := i\v((a — a}) (cr+ D(irj) + er_ D{irj)' { ) — irj (o + D(ifj) — er_ D(irtfy) 
= i\v(a-a { )(<j + + <7-) (21) 

+ * A ^E ^ T -(a 2 -a t2 + l-m)(a + at)™-i( (T+ + (_i)™ (J _). 

m— 1 
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Notice that the 'balanced detuning' 5 is, unlike 6, always positive so that the 
degeneracy condition for Hq is now: v = mS, m = 1,2, . . . . 
The unitary operator Ta has the following explicit form: 



T \ 



xX D{ie A rj) k a D{ie\ rf) 
-x A D(ie+T 1 Y x+Diie^rjY 

x+D[i{r,-n)/2) x A D{i{f, + rj)/2) 
-x A D(i^ + r,)/2)^ x+D {i(r, - V )/2) ] 



(22) 



with 



1 



A 



'1 



4 2V4 + A 5 



1 



± sign(A) . - 



1 



4 2V4 + A 2 



1 

± - 



(23) 



(24) 



A 2V4 + A 2 2 

The unitary operator Ta can be decomposed as the product of three unitary 
transformations: 

T A = T 3 T 2 T 1 . 

The transformation T\ has been introduced by Moya-Cessa et al. [35] and has 
the following form: 



\/2 V2 



1 

72 



-2?t £> 



D = D(iri/2). 



The transformation T 2 is a spin rotation by the angle round the y-axis: 



T 2 := 



cos 0/2 -sin 0/2 
sin 0/2 cos 0/2 



where the angle 0, — 7r/2 < < 7r/2, verifies the condition 

2 2fi fl 
The third transformation is given by 

T3 := \ (D(irj/2) + D(im f ) + \ {Diiv/2) - D{if,/2)^) a z 

D{ifi/2) 
D(irj/2)i 
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One can check that in the weak field limit SI;? — > (A = 6/CIr ~> ±00) and 
in the strong field limit — * +00 (A — * 0) the transformation Ta has the 
following behaviour: 



lim Ta = Id, lim Ta — ?2(— 7r) = 

A— >+oo A^ — 00 





-1 



(weak field limit); 



lim Ta 

A^O 



(strong field limit). 



Thus, in the weak field limit, Ta goes to the identity if S > and to a spin 
rotation which sends a z into — a z for S < (so that, as already observed, 5 is 
always positive); while, in the strong field limit, it goes to the transformation 
introduced by Moya-Cessa et alii. 

The constant A which appears in the expression of the interaction component 
Hf of the balanced Hamiltonian H will play the role of perturbative parameter 
in our later analysis. Moreover, it will be seen that the constant 77 rules the 
number of terms that have to be considered at each perturbative order. We 
remark that both A and fj are bounded functions of the Rabi frequency £Ir; 
indeed: 

0<1A]= = M<~M 



6 2 



and 



\s\ 



\v 



< 



>1\ 



Thus, one can apply a perturbative approach also in the case of a large Rabi 
frequency (hence in presence of an intense laser field). 



4 Perturbative analysis: outline of the method 

Let Sjo, 55 j be hermitian operators and assume that 9)q has a purely discrete 
spectrum (i.e. it has a pure point spectrum with finite-dimensional eigenspaces) . 
Denote by 

Eq < Ei < E2 < . . . 

the (possibly degenerate) eigenvalues of $) and by Po, Pi, P2, ■ ■ ■ the associated 

eigenprojectors. 

Now, consider the operator 

£(A)=£ + A% AeC, 

which is hermitian if A is real. It is possible to show that, under certain con- 
ditions (see |37| |38p. there exist positive constants ro,ri,T2, • ■ ■ and a simply 
connected neighbourhood X of zero in C such that the following contour integral 
on the complex plane 

P m {\) = ^-ldE (ftW-Ey 1 A eX, (25) 

^ J\E-E m \=r m 
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defines a projection (P m (X) 2 — P m (A)), which is an orthogonal projection for 
real A, with P m (0) = P m , and 19 An P m (A) is an analytic operator- valued 
function. Moreover, the range of P m (A) is an invariant subspace for f)(A), hence 

$)(\) P m {\) = P m (A) 55(A) P ro (A), (26) 
and there exists an analytic family U (A) of invertible operators such that 

P m = U(\)P m (X)U(X)-\ U(0) = ld, (27) 

and 

U(X) = e lZ ^ A el, (28) 

with Z(X*) = Z(Xy (hence, for real A, Z{X) is hermitian and U{X) is unitary), 
where I 3 X i— > Z(X) is analytic. Let us observe explicitly that the following 
relation holds: 

U{X*) = e jZ(A * } = e lZ ^ = U(X)-^. (29) 

We remark that, as it is easily shown, the function A <— > U(X) is not defined 
uniquely by condition l|27|) even in the simplest case when f)o has a nondegen- 
erate spectrum. 

Now, let us define the operator £j(A) by 

^(X):=U(X)^(X)U(X)-\ (30) 

which, for real A, is unitarily equivalent to -f)(A). Using relations (|26|l and (|27(l . 
we find 

5(A) P m - U(X)^{X)P m {X)U(X)- 1 

= C/(A)P m (A)jo(A)P m (A)[/(A)- 1 

and hence: 

5(A) P m = P m 5(A) P m m = 0, 1, . . . . (31) 

It follows that 

[fa, 5(A)] = (32) 
and then we obtain the following important decomposition formula 

U(X) £(A) ^(A)" 1 =f)o + C(X), (33) 

where [C(A),.f)o] = 0, i.e. C(A) is a constant of the motion for the evolution 
generated by Sjo. Notice that by virtue of relation Ij29(l we have: 

C(A*) = C/(A*)^(A*)J7(A*)- 1 -fx, 
= C/(A)- 1 t^(A) t C/(A)t-^ 

= C(X)K 

Thus, in particular, for real A, C(A) is hermitian. 
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At this point, we are ready to build the perturbative decomposition of the 
operator f)(A) anticipated in section ^ Indeed, since the functions A ^ C(A) 
and A Z{\) are analytic in I and C(0) = Z(0) = 0, we can write: 

oo oo 

C(A) = A m C m , Z{\) = J2 A m Z m Ael; (34) 

m— 1 m— 1 

here we notice that — as C(A*) = C(A) 1 * and Z(A*) = Z(X)' { — for any m, 
C m = C m and Z m = Z m . Then, if we set 

4 n) (A) := e -i(A2i+-+A»z.) aoe i(XZ 1 +...+A%) i (35) 

£ ( n>(A) . = e -i(AZ 1+ ...+A»Z„)( ACi+ .,. + r g e i(AZ 1 +...+A»Z B ) ) ( 3 g) 

by virtue of formula (|33(l we find that equation is satisfied together with 
condition J2J. 

Now, in order to determine the operators {C„(A)} and {Z n (X)}, let us substitute 
the exponential form e tZ ^ of t/(A) in formula (|33|l : we obtain: 

+ E 5 ad ^ ^ (A) = ^° + C(A) ' (37) 

where we recall that adx L = [K, L]. Next, substituting the power expan- 
sions (|34f) in this equation, in correspondence to the various orders in the pe- 
rurbative parameter A, we get the following set of conditions: 

Ci-i[Zi > fl ]-%=0 , [C 1 ,Sj }=0 (38) 
C2-i[Z 2 ,S)o] + l[Z 1 ,[Z 1 ,Si ]]-i[Z 1 ,Si t ]=0 , [C 2 ,i3o] = (39) 



where we have taken into account also the additional constraint [C(A),^o] = 0. 
The generic term, after the first one, in this infinite sequence of equations is 
easily shown to be the following: 

°n ~ E E [Z kl ,[..-,[Zk m ,Sj }...}} 

^— ' ml z — ' 

m=l feiH \-k m =n 

- E S E [^ i; [...,[^ m ,^]...]] =o, 

m=l ' k\-i \-k m =n— 1 

[C„,i5o] =0 n = 2,3,... . 

The infinite set of equations can be solved recursively and the solution, as al- 
ready anticipated, is not unique. The first equation, together with the first con- 
straint, determines Z\ up to an operator commuting with and C\ uniquely. 
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Indeed, since 

oo 

C x = P m Ci P m and [Zi , So] = 2 (E k - Ej ) Pj Z x P k , (40) 

m=0 j=jLk 

we conclude that 

oo 

C 1 =^ l P m Si l P m (41) 

m=0 

and 

oo 

Zi = ^ P m Zi P m + t ^ (£7 fc - P J % ft- ( 42 ) 

m=0 j/fe 

This last equation admits a 'minimal solution' which is obtained by imposing a 
further condition, namely 

P m Z 1 P m =0 m = 0,l,... . 

For n > 1, we will use an analogous reasoning. Indeed, let us define, for n > 2, 
the following operator function: 

" i 7 " 



G n (Zi, . . . , Z„_i) :— — - [Z fel , [. . . , [Z km ,Sj ] . . .]] 

m=2 fexH hK m =n 

-m 

+ E ^ E [^,[...,[^, %]•■•]]• (43) 

m— 1 fciH hfe m = 71— 1 

Now, assume that the first n equations have been solved. Then, the operator 
G n+ i(Zi, . . . , Z n ) is known explicitly and hence 

oo 

= ^ m ■ ■ ■ j -^n) ^m) (44) 

m— 

[Zn+i,So] =i^P J -G n+ i(Z 1 ,...,Z n )P fc . (45) 

Again, this last equation determines Z n+ \ up to an operator commuting with 
fio- In general, the choice of a particular solution for Z n+ \ will also influence 
the form of C n+2 , Z n+2 , .... 

Thus, we conclude that the sequence of equations defined above admits infinite 
solutions (even in the case when S^q has a nondcgenerate spectrum). This fact 
had to be expected as a consequence of the non-unicity of U(X). Anyway, there 
is a unique 'minimal solution' {Ci, Zi, . . .} which fulfills the following additional 
condition: 

P rn Z k P m = m = 0, 1, . . . , k = 1, 2, . . . . (46) 
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The preceding scheme can be extended to the general case when the inter- 
action component Joj(A) of Sj(X) does not depend linearly on A: 

oo 

£(A) =£o+%(A) = 5> m .f) m . (47) 

Indeed, given an operator X, let us set Fo(X) = X and 

F n (X;Z 1 ,...,Z n ) := V — V [. . . , [Z km ,X] . . .]] , (48) 

A — ' ml z — ' 

m=\ k\-\ \-k m =n 

for n > 1. Then we can define the operator function 

n 

G n ($Jo, . . . ,Sj n ; Z\, . . . , Z n —i) := S ' F n ^ m (Sj m ; Z%, . . . , Z n _ m ) 

m=0 

- *[2„,i3o] »>1, (49) 

which generalizes definition (|43|) . At this point, one can show that this time the 
decomposition formula (|33[> leads to the following sequence of equations: 

C n — i[Z n ,f)o] = G n ($jo, . . . ,f) n ; Z%, . . . , Z n ^i), 

[C n ,S) ] = neN. (50) 

Again, the general solution of this set of equations can be obtained recursively 
by formulae (|44|l and H45|) . We want to show next that it is possible to give 
simple explicit expressions for the operators C ra , Z n which do not involve the 
eigenprojectors of fjo- 

To this aim, let us denote by Ho, Hi,... the eigenspaces associated with the 
eigenvalues Eq, E%, . . .; namely, let us set: 

H n = Ran(P n ), d n = dim(7i„) n = 0, 1, . . . , d_i = 0. 

Now, let {\n) : n = 0, 1, . . .} be an orthonormal basis formed by eigenvectors 
of i^o such that 

71-1 

H n = Span(|c„), . . . , |c„ + d n )}, c„ = ^ d n . 

fe=-i 

We will denote by A, , N respectively the annihilation, creation and number 
operators relative to this basis: 

A f \n) = Vn + 1 \n + 1) , N\n) = n\n). 

Then we can write 

f>o = Efi (51) 
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and express the operators G\ = G n = G n ($)o, . . . ,S) n ; Z\, . . . , Z n -i), n > 2, 
in the following form: 

OC 

G n = (N) + (9 l n ] (N) A m + At >»gM (N)) . (52) 

ra— 1 

At this point one can check that the operators 

d-l 

Cn = g [ ?J (AO + (x ( E N +m ~ E n) 9 l n ] (N) A m + h.c.) , (53) 

rn — 1 

oo 

Z n = Z® + iJ2 (7 (E$ +m Efi) (N) A m - h.c.) (54) 

m= 1 

- where (i = sup{ d„ : n = 0, 1, . . .}, Z$ is a hermitian operator such that 
[z[?',flo] = and the functions Xil are defined by 

X(0) = 1, *(*) =0 a; ^ 0, 

7 (0) = 0, -fix) = - x^O 
x 

— are solutions of eq. (| 5 C)|> . This can be readily verified by direct substitution 
if one observes that 

X(x)x = 0, j(x)x = l-x(x). 

We remark that if now {\\n} : n = 0, 1, . . .} is any orthonormal basis of eigen- 
vectors of S)o, so that 

9) \\n) = E(n) \\n) , 

where in general E(n) 7^ E n , and A, A* , N are respectively the annihilation, 
creation, and number operator associated with this basis, then formulae l (5*5j l 
and (|54|l still apply with the substitutions 

E^ 1 — ► E(N), d-l 1 — ► 00. 

We observe explicitly that with these modifications formulae (|53f) and (|54|l can 
be applied formally ignoring the possible degeneracies in the spectrum of .fio- 

5 Perturbative analysis: treatment of the BH 

In this section, we want to apply the perturbative method outlined in the preced- 
ing one to the balanced Hamiltonian H . More precisely, for better illustrating 
the method, we will use, beside H, a Hamiltonian which is unitarily (hence 
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physically) equivalent to H . 

Indeed, let us consider the unitary operator 



/ 7T 

T := exp\i-na a 



"K 

cos I — n 



sin (~n)o- x - (55) 



.2 

One checks easily that the following transformation formulae hold: 

TaT T = — iaa x , (56) 
Ta z T^ = j™a z , (57) 

Ta±T^ = cos 2 (|n) cr± + sin 2 (| n) <r T . (58) 

With the aid of this formulae, we find that 

H :=THT* = H + H t , (59) 



where: 



Ho := vh + ^Se inA a z + Xfju (60) 



and 

H t := \is(a + a r ) (61) 

+ \ v V" TL ( fl ta _ fl 2 + ! _ m \ ( fl t _ a jm-l g *m™ , + (_1)«V + ) . 



m— 1 



Now, if we assume that 77 <C 1, we can drastically reduce the number of terms 
that need to be considered in the interaction components of the Hamiltonians 
H and H: namely, we can set: 

H x ~ i\v (a- a 1 ) (er + + er_) + A77Z/ (a t2 - a 2 ) (cr + - cr_) (62) 
#j ~ At/(a + a t )+Ar;i/(a t2 -a 2 )e l7rft (CT_-(7 + ). (63) 

Moreover, in order to treat H-t and H-t as perturbations, with perturbative 
parameter A, we will also assume that A <C 1. For simplicity, we will set 

H = e -<^ + ^+-) (H + C 1 + C 2 + ---)e i (^+^+-), (64) 

H = e -<^ + ^ + -) (H + C 1 +C 2 + ---)e<^ + ^+-). (65) 

Thus, differently from section the perturbative parameter A will be included 
in the operators C n , C n , Z n , Z n . 

Our purpose is to perform the perturbative analysis up to the second order; 
precisely, since we want to compute the first and second order corrections to 
the energy spectrum and the first order expression of the evolution operator, we 
need to compute the operators C\, Zi, and C%. 
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At this point, it is convenient to distinguish three cases: when r\ is much smaller 
than A, when rj is of the same order of magnitude of A, when 77 is much larger 
than A. 

Let us consider the case when r)< A. 
In this case, the second term which appears respectively in the r.h.s. of eq. I|62[l 
and eq. , due to the presence of the factor A 77 <C A 2 , can be skipped in a 
second order treatment; hence we can set 

if j ~ i A v {a — a^) (er + + a J) , i?j ~ A v (a + a 1 ) . 

Let us use the Hamiltonian H first. Observe that, within the given approxima- 
tion, the subspaces H g , H e of the total Hilbert space H, H g © Ti e = H, 

H g := Span{|n) ® \g) : n = 0, 1, . . .}, H e := Span{|n) ® \e) : n = 0, 1, . . .}, 

are invariant subspaces for H. Thus, it will be convenient to split the perturba- 
tive problem by restricting to each invariant subspace. In order to do this, let 
us denote by a g , n g , H g , Hq 7 Hi (resp. a e , . . .) the restriction of the operators 

a, n, H, Hq, i?| to 7i g (resp. to 7i e ). Then, we have: 

H g = H° +Hf = vh g - l -5 e-" S + A v (a g + a\) . 

At this point, identify the basis {\n) ® \g) : n = 0, 1, . . .} in TL g with the basis 
{\\n) : n = 0, 1, . . .} of section 0] and the operators a g , n g with A, N. Then, 
applying formula we find: 

Cf = A v X (y + 5 e l7r ™ 9 )a g + h.c. . 

Hence, Cf — unless the resonance condition v — 5 is satisfied; in this case, we 
have: 

ei = \\v{a g +a\) + \\v{a g -a\)j™°. 

In the subspace H. e , arguing as above, we find that Cf is zero unless v = 5, in 
which case: 

C\ = \\ v (a e + at) - l -\ v (a e - oj) e 1 ^ . 

Hence, considering the total Hilbert space TL, at the first perturbative order in 
A we find that 

C x = for u^5, (66) 
C x = ^Ai/(a + a f ) --Xv(a-a i )e i ' !Tn a z forv = 5. (67) 

Moreover, one can check that the minimal solution for Z\ is given by: 

Zt = -\v~f{v-8) {{a -a)) ~ (a + a)) e l ™ a z ) 

+ l -\^((a-J) + ( a + J)e l ™o- z ). (68) 
^ v + 8 
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The operator Z\ allows us to compute the second order constant of the motion 
C-2- The result is found to be: 



Co 



v + 8 W 2 / 2 
- X^ 2 7 (u-S)((n+^je i ^a z + ^). (.!!)) 

On the other hand, in the 'reference frame' associated with H , at the first order 
we have: 

Ci = T t (7iT = for v^S, (70) 
C\ = iXu — a <t_) for u = 8. (71) 

Notice that, in correspondence to the resonance condition v = 8, C\ coin- 
cides with the constant of the motion S, with = 5, of the Hamiltonian 
Hjc which is unitarily equivalent to the Jaynes-Cummings Hamiltonian (see 
section |2J) and hence with the prescription of the RWA (see section |3J) . We 
could have obtained this result using directly the Hamiltonian H, observing 
that mv = 8, m = 1,2, . . ., is the degeneracy condition for Hq and then apply- 
ing formula (|41|l . 

Thus, it could seem that, in the resonant regime v = 5, the first order approxi- 
mation of H coincides with the RWA Hamiltonian Hq + C\ . Actually, it would 
be so if Z\ was identically zero (recall formula (|64|) ~). but we have: 

Z 1 = T^Z X T 

= — \v^{y — (5) (aa + + a) er_) — A (a er_ + a 1 er + ) . (72) 

z/ + 5 

We will see the effects of this fact in section [7| 
At the second order, we find: 

C 2 = T^C 2 T 

V 1 ' ■ hi) - - I) - - J) ( (» + i) c. + i). (73) 

Let us now consider the case when 77 ~ A. 
In this case, since A 77 ~ A 2 , we have to consider the second term in the r.h.s. 
respectively of eq. (|fi2"|) and eq. (SHjl at the second perturbative order. In this 
case, one can work directly with the Hamiltonian H. Then, the result is the 
following. At the first order, the constant of the motion C\ is still given by 
formulae (|70|l and l|71|) . Z\ by formula l|72|) . At the second order, one finds: 

4 - A2 TTi ((:• + 5) -5) - A2 " 27( " - 11 (( ft + i) " + i 

- A?7i/x(2^- 8) (a 2 a + + a f 2 a J) . (74) 
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Hence an extra term, with respect to formula (|73|) , appears in the expression of 
C%. This term is related to the resonance 2v = 5. 

Eventually, let us consider the case when rj 3> A. 
In this case, we have: 

A 2 « Aij « A < ij « 1. 

Hence, the term proportional to A fj can now be regarded as the leading term 
at the second perturbative order. Then, the operators C±, Z\ and C% are given 
again by formulae l(7Djl. J7TJ, JUJl and (f7^|l . 

A special attention is deserved by the nearly resonant regime. Indeed, if the 
condition \v — S\ <C v is satisfied, it is convenient to set H = Hq + H'^, where: 

H' Q = v (n + ^a^j +\rjv, (75) 

+ i \v [a - a f ) (o+ + aJ) + Xfjv (a f2 - a 2 ) (cr+ - cr_) . (76) 

Thus, the zeroth order Hamiltonian iJp has degeneracies in its spectrum as in 
the resonant case and an extra term appears in the perturbation. This strategy, 
which is analogous to carefully choosing the origin of the power expansion of an 
analytic function, allows to obtain larger convergence radii for our pertubative 
expansions. In the nearly resonant case, one finds that the expressions of the 
operators C\, Z\ and (% are the following: 

C% = - (5 — v) a z + i A v (a cr+ — a? <r_) , (77) 



Zx = -- A (a<7_ + a f o- + ) , (78) 




Notice that, as in the exactly resonant regime, the operator Hq + C\ coincides 
with the RWA Hamiltonian (now obtained using as reference Hamiltonian Hq), 
but again, since Z\ 7^ 0, it is not the correct first order approximation of H. 



6 Corrections to eigenvalues: the RWA and the 
Bloch-Siegert shift 

By means of the perturbative method discussed in the previous sections one can 
obtain approximate expressions, at each perturbarive order, of the eigenvalues, 
the eigenprojectors and of the evolution operator associated with the BH. 
For instance, recalling formulae l|77() and l|79l) . one finds that approximate ex- 
pressions, at the second perturbative order, of the energy levels of the BH in 
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the nearly resonant regime, i.e. for \v — 5\ -C v 1 are given by the eigenvalues of 
the hermitian operator 



I) 



(2) _ 









(n-\ 











+ i A v (a a+ - a) cr_) + i A 2 i^ f f n + er 2 - i j . 

Observe that the eigenspaces of the operator f)( 2 ) are the one-dimensional sub- 
space Hq = Span{|0) ® |<?)} and the two-dimensional subspaces 

U n = Span{|n- 1) <g> |e), |n) <g> n = l,2,.... 



Hence, all one has to do is to diagonalize the matrix representation of \) 2 ' in each 

of these two-dimensional subspaces; namely, one has to diagonalize the 2x2 

matrices _ 
2l n + <B n i A v yfn 

-iXv^/n 2l„ - 25„ 

where: 



ra = l,2,. 



(80) 



»„ = 



+ Ary^--A z/, 



A 2 ^ n. 



(81) 
(82) 



Thus, the energy levels of the BH H, in the nearly resonant regime, are given 
by 

1 1 2 

i?o — — f + A fj v A v, 



E n ±~v[n--j+\r } v±^-((5-v) + \ 2 vn S j +X 2 v 2 n--\ 2 v : (83) 

with n = 1, 2, . . . . The previous expressions coincide, skipping the second order 
corrections, with the ones obtained applying the RWA. This is due to the fact 
that, as already observed in section the operator Hq + i A v {a a + + aV_] 
coincides with the result of the application of the RWA to the Hamiltonian H, 
in the nearly resonant regime \v — 5\ <C v. Hence, on the one hand, the RWA 
gives the correct first order expressions for the eigenvalues. In the next section 
we will show that, on the other hand, the RWA does not give the correct first 
order approximation of the evolution operator associated with the BH (hence, 
with the ITH). 

There is also a second order effect that cannot be predicted if one simply 
applies the RWA. This effect may be compared with the so called 'Bloch-Siegert 
shift'. In order to clarify this point, let us consider the matrix representation of a 
hermitian operator () (which can be thought as a perturbative approximation of 
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the Hamiltonian of a physical system) in a two-dimensional invariant subspace, 
with respect to an orthonormal basis {|1), |2)} in this space: 



M = 

Now, if one sets 



f)n f)l2 
f)21 f)22 



f)H = f)22 = f)22. f)l2 = f)21- 



21 = 



1)11 + f)2 



b = 



- f)22 



C = f)!2 = f>21) 



the matrix M can be rewritten as 



M 



21+ b c 
c* 21- b 



and the eigenvalues of the hermitian matrix M are given by the following simple 
formula: 

E±(b) =2l± ^/b 2 + \c\ 2 . 

Let us assume that t ^ 0. Then the vectors |1),|2) are not eigenvectors of 
f) and the graphics of the functions b >— > E±(b) are the two branches of a 
hyperbola whose asymptotes intersect at the point of coordinates (0,21). Thus 
the difference between the two eigenvalues E + (b) — E-(b) attains its minimum 
('level anticrossing') at b = 0. This is also the condition for which the transition 
probability Pi_>2(£) assumes periodically the value 1 (otherwise Pi^it) < 1). 
Indeed, according to a well known formula, we have: 



P^ 2 (t) := |(2|ex P Hf)i)|l)| : 



Vb 2 + \c\ 2 t 



At this point, suppose that — due, for instance, to (higher order) perturbative 
corrections — the matrix M undergoes a modification of the type 



M 



M 



Then the level anticrossing condition undergoes a shift: b + e = 0. In our case, 
we can do the following identifications: 



2l = 2l„, b = ^[d 



~ \ 2 vn, 



1,2,. 



Thus, the second order level anticrossing shift for the subspace Ti. n is given by 

Mr 



1 ^2 1 
- A vn = - 

2 



n= 1,2,. 



2 y^TP 

A similar phenomenon appears in the classical work of Bloch and Siegert [3T] 
on the magnetic resonance, whose Hamiltonian can be replaced, using Floquet's 
theorem, by a time-independent effective Hamiltonian containing virtual terms, 
as shown later by Shirley [35]. The presence of these virtual terms gives rise, 
at the second perturbative order, to a level anticrossing shift, which translates 
into a shift of the magnetic resonance condition, the Bloch-Siegert shift. 
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7 The evolution operator 

From this point onwards, for the sake of conciseness, we will use the following 
notation. Given a couple of functions / and h of the perurbative parameter A, 
if /(A) = h{\) + O^X 2 ), we wil1 write simply 

/(A) £ h{\). 

Then, let il(i, t ) be the evolution operator associated with the ion trap Hamil- 
tonian H(t): 



d 

i— ii ) (t, t ) = H(t) U(t, t ), il(*o, to) = Id. 



As we have seen, il(t, to) can be decomposed as 

il(t, t ) = R\ Tl e -* A (*-*o) T A . 

Moreover, the evolution operator associated with H, namely (£(i) 
admits a perturbative decomposition. At the first order, we have: 



(84) 
(85) 

-iHt 



<£(t) ~ cxp (-t e"^ 1 (# + #i) e^ 1 t) 

_ e -iZi e -i(H +C!)t e »Zi 



= e 



-iZi 



e lZl =: <£i(*), 



(86) 



where for obtaining the third line we have used the fact that [C\,Ho] = 0. In 
order to deal with simpler formulae, let us consider the case when the resonance 
condition v = 5 is exactly satisfied. In this case, we have that Ho = v (n + ^ <r z ) 
and C\ = i\v (clcf + — a* a-). Thus, C\ is, up to a unitary transformation, the 
infinitesimal generator of the Jaynes-Cummings evolutor JO(t); hence: 



JC(t) := exp l-idt 



^ A JJ(t) e** ft 
cos (A v \Jh + 1 t) 



sin (A v \Jn + 1 t) 



Vn+1 



sin ( A i> i ) / \ 
— — — at cos ( A v t ) 



(87) 



Then it turns out that the expression of the evolution operator 9i(t) obtained 
applying the RWA to H, namely 



«R(i) = exp(-iH t) JO(t), 



(88) 
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would be correct at the first perturbative order only if the operator Z\ was 
identically zero. But this is not the case since, for v = 8, it turns out that 
Z\ — — i A (a <t_ + a) a + ) and we have: 



exp 



cos I I A \fh 



sin 



sin(iAVn+l) n ~r\ 

-i z -a cos(5AVn-+l) 



\/n + 1 



This result can be also expressed saying that the RWA neglects the first order 
correction to the unperturbed eigenprojectors. 

In order to get a more explicit comparison of the correct first order ap- 
proximate evolution operator with the RWA evolution operator £H(i) , we 
proceed as follows. First, we observe that 



d(t) 



e -iH t jj^ e iZ t 



-iZi -iH t iZ 



1 Jb(t). 



Then, using formula (|89|l . we find that the operator e 
following expression: 



■iZ\ ^ — iHnt c iZ\ 



has the 



(a x {h) + p x {n)e i2vt ) e^'M)* Kx (fi) (1 - e l2vt ) e^^+i)*^ 

K X (n + l) (1 - e~ avt ) e'^-i^a (a A (n+l) + /3 A (n + l)e- 42l/t ) e^K^*)* 



where we have set 



a\(n) :— cos 2 ( — A \ffi I , (3\(n) ■— sin 2 ( — A \ffi I , 



cos^i A Vftj sin^i A y/foj 



Thus, due to the fact that K\(h) 



■ Aid, we have: 



g- <z i e -iH t & iZ 1 *t e -iH„t 



for i ^ 0. 



It follows that, for t ^ 0, 



This proves that the RWA does not provide, already at the first perturbative 
order, the correct approximate expression of the evolution operator associated 
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with the BH (hence with the ITH). 

Now, in order to obtain a direct comparison of €(t) with 9t(i), observe that 
<£(t) ~ e^ 1 e~ lflot e^ 1 JJ{t) 

= e-^W e^-V e-*™ JJ(t) 

£ e -i(Zi(0)-Z 1 (-t)) e -iffot i £f( t ) j 

where t <—> Z± (t) is solution of the Heisenberg equation 

(JfcZ 1 > \(t) = -i[Zi(t),H ], Zi(0) = Zi. 

Then, since 

-i[2i(-t),#o] = iAi/(a<7_e i2vt + a t <7 + e- i2, ' t ) = : Yi(f), 
we find the following formula: 

<£(*) ~ exp^-i^ Yi(r) dr^ SR(i). (90) 

This expression provides a direct relation between <£(i) and $H(t). Notice that 
it contains the integral of an oscillating function. Anyway, since this integral 
appears as the argument of an exponential, we are not led to the erroneous 
conclusion that its contribution can be neglected, as it is often incorrectly argued 
using a Feynman-Dyson expansion of the interaction picture evolution operator 
(see section EJl. 



8 Discussion 

In writing the present paper, the authors had in mind two main aims: 

• to show that the Hamiltonian of a trapped ion interacting with a laser 
field can be studied even if the condition 17^ <C v is not satisfied, since 
this condition is incompatible with applications that physicists consider 
to be relevant nowadays, for instance fast ion trap quantum computers; 

• to show that a rigorous perturbative approach can improve the results 
obtained by simply applying the RWA, still preserving the chance of per- 
forming explicit and manageable calculations. 

With regard to the first point, it has been shown that the study of the ion 
trap Hamiltonian can be reduced to the study of a time-independent effective 
Hamiltonian (the BH) in which the coupling constant is scarcely sensitive to the 
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Rabi frequency Qr, hence perfectly fit for our aims. The resonance condition 
for the BH has a very simple form: 



m v = 5 = 



AVt 2 R + 5 2 , m = 1,2,... , 



where we recall that S denotes the ion-laser detuning. 

With regard to the second point, we have shown that a suitable perturbative 
approach allows to write a very powerful perturbative expansion of the evolution 
operator of the system. Indeed, using the notation of section 0] we have: 



where the operators Ci, Z%, Ca, Z%, . . . can obtained by a recursive algebraic pro- 
cedure. Notice that any truncation of this perturbative expansion is a unitary 
operator, a very valuable feature for an approximate expression of the evolu- 
tion operator. Our approach, whose general validity goes beyond the argument 
of this paper, can be applied successfully to the BH pointing out two main 
facts. First, the RWA terms and the counter rotating terms play different roles 
(but with the same dignity) in the perturbative expansion (|91|l . The former 
appear in the time-dependent component of the expansion (the one associated 
with the operators Ci, C2, . . .), while the latter appear in the time-independent 
component (associated with the operators Z\, Z2, ■ ■ .). Second, already at the 
first perturbative order the counter rotating terms, completely neglected by the 
RWA, give rise to a correction that can be regarded as a perturbative correc- 
tion to the unperturbed eigenprojectors (while, as we have seen, the first order 
correction to the eigenvalues coincides with the prescription of the RWA). In 
conclusion, we believe that our approach can provide more accurate expres- 
sions of the evolution operator of the ion trap Hamiltonian for a wide range of 
intensities of the driving laser field. 
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